Two freeze-out model for the hadrons produced in the Relativistic Heavy-Ion 
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An expanding fireball model with two freeze-outs, which assumes that the chemical freeze-out 
occurs earlier at higher temperature and the thermal freeze-out occurs later at lower temperature, 
is developed and successfully applied to fit simultaneously the hadron ratios, transverse momentum 
spectra of measured hadrons and rapidity distribution of charged hadrons measured in Au-|-Au 
collision at 200 A-GeV by STAR, PHENIX and PHOBOS collaborations. The quality of all the 
fittings is very good and the resulting parameters are in agreement with published values. 

PACS numbers: 24.10.Pa,25.75.-q 



I. INTRODUCTION 



The hadron production in relativistic heavy-ion colli- 
sions have been extensively studied: The multiplicities of 
various hadrons or the ratios among them are nicely fitted 
using the statistical distributions with a few parameters 
such as the temperature, Teh, the baryon chemical po- 
tential, fiB.ch, the strangeness chemical potential, ^s,ch 
and the strangeness fugacity, 7s[l-6]. The temperature 
thus obtained at RHIC energies is so high and close to 
the phase transition temperature to quark-gluon plasma 
and suggests that the hadrons are chemically frozen out 
just after the hadronization. The value of 7^ close to 
one suggests that the strangeness is nearly equilibrated. 
The slopes of the transverse momentum spectra below 
2GeV/c of various hadrons are also nicely fitted from an 
expanding fireball model with a single set of parameters 
such as the temperature Tth, the baryon chemical poten- 
tial ^B.th, the strangeness chemical potential fis,t!n and 
the transverse expansion velocity /3[7]. The success of 
the thermal analysis with the Lorentz boosted thermal 
distribution is regarded as the evidence of the radial ex- 
pansion. 

However, the temperatures of the two analysis are dif- 
ferent, i.e. Teh 7^ Tth and thus one cannot fit both 
the magnitudes and slopes of the transverse momentum 
spectra of various hadrons. In other words, when one 
fits the transverse momentum spectra normalization con- 
stant has to to be adjusted separately for each hadron 
species. The difference in the two freeze-out tempera- 
tures are interpreted as the chemical freeze-out occurs 
earlier at high temperature where the inelastic collisions 
of the type A-l-B ^ C-l-D becomes less frequent and the 
numbers of each hadron species are no more changing and 
thus kept fixed, while during the continuing expansion 
of the system, elastic collisions among the same species 
maintains the thermal equilibrium and lowers the tem- 
perature until the thermal freeze-out. 
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In the hydrodynamic calculations, the difference of the 
two freeze-outs has to be taken into consideration when 
one wants to compare the results of calculation with the 
experimental data. Such models reported so far are the 
hybrid models[9-ll] and the partial chemical equilibrium 
model(PCE) [11-15]. In the hybrid models or the hy- 
dro-ftransoport models, the expansion is described by 
the hydrodynamics, starting from an initial state of QGP 
and at Tc the equation of state is changed to the hadronic 
EOS. At a certain switching temperature in the hadronic 
phase, Ts^, hadrons are generated by the Monte Carlo 
method and the subsequent interaction among hadrons 
is described using microscopic hadronic transport models 
such as UrQMD[9, 10] or JAM[11]. In the PCE model, af- 
ter the chemical freeze-out the numbers of hadrons whose 
lifetime is longer than the characteristic time(~ 10 fm/c) 
are kept fixed, while the reactions involving hadrons with 
shorter lifetime(e.g. A ^ p + tt) are still in chemical 
equilibrium. In this case one has to solve large number 
of coupled equations to get the values of chemical poten- 
tials. 

In the blast-wave model the radial expansion is one 
of the essential ingredients in fitting the transverse mo- 
mentum spectra with the Lorentz boosted thermal dis- 
tribution which does not change the total number. Thus 
in principle, if one uses the Lorentz transformed distri- 
bution functions in the chemical analysis, the result is 
expected to be the same as that without the expansion. 
Then it seems natural to make an expanding fireball 
model in which the chemical freeze-out occurs earlier at 
higher temperature and then thermal freeze-out occurs 
at lower temperature, while in between the two temper- 
atures the number of each hadron species is fixed. Even 
though the concept of the two freeze-outs are already 
adopted in the hydrodynamic calculations as discussed 
above, there is no blast-wave type calculations reported 
until now which fits both the ratios and the particle spec- 
tra of various hadrons in a single model. Then the range 
of parameters which fit the data with similar quality in 
the chemical and thermal analysis may be pinned down 
to a single set of parameters. 

In this work we report on a simultaneous blast wave- 
type analysis of both the hadron ratios and transverse 



momcntuin spectra measured at RHIC assuming that the 
chemical freeze-out occurs earher at higher temperature 
and the number of each hadron species is kept fixed un- 
til the thermal freeze-out when the particles escape from 
the system to reach the detectors. At the temperature 
higher than the chemical freeze-out temperature Teh, the 
chemical potentials for hadrons are the algebraic sums 
of the baryon and strangeness chemical potentials. The 
ratios among various hadrons are fitted using the blast- 
wave equation for the transverse momentum spectra [16] 
by integrating over pt. The contribution from the de- 
cay of high lying resonances are also taken into account 
by integrating over pt the resonance contribution to the 
transverse momentum spectra[18]. 

Below Teh, the chemical potential of each hadron 
species is no more the linear sum of the baryon and the 
strangeness chemical potentials but it should be numer- 
ically found to guarantee the fixed thermal number of 
hadrons. Using the chemical potential for each hadron 
species thus obtained one can fit the momentum spectra 
as done by Dobler et al[16] which is the most sophisti- 
cated version of blast-wave model. 

In this way we have successfully fitted both the ra- 
tios among hadrons, transverse momentum spectra of pi- 
ous, kaons and protons, and the rapidity spectra of the 
charged hadrons. 

One difference of this calculation with the PCE model 
is in the treatment of the short lived resonances. In this 
work they are treated as stable particles until the ther- 
mal freeze-out and their decay has been considered after 
the freeze-out. This approximation causes a small errors 
but reduces the computation significantly since the cou- 
pled equations for the chemical potentials now reduces 
to independent equations. 



II. BLAST- WAVE MODEL WITH TWO 
FREEZE-OUTS 



rmaxiv) = -Roa/I 



but in this work it is kept con- 



stant. Then one gets for the hadron spectrum 
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where 7s is the strangeness fugacity. 

It should be emphasized that the chemical potential, 
pn is different above and below the chemical freeze-out 
temperature, Teh- Until the chemical freeze-out (for 
T > Teh), Pi is the algebraic sum of the baryon and 
the strangeness chemical potentials: 



T 



Pi = (nq - nq)pq + {us - ns)ps 



(4) 



Below Teh, however, Eq. (4) does not hold and pi should 
be obtained from the number of thermal hadrons of type 
i, which is fixed at the chemical freeze-out as will be 
discussed later. 



For the number of thermal hadrons, 
grates Eq. (3) over mTdrriTdy. 
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The total number of hadrons measured is the sum of the 
thermal ones and the decayed ones from other hadrons 
or resonances. The resonance contribution to the hadron 
momentum spectra is nicely treated by Sollfrank[18], and 
to get the number of hadrons from the resonance decay 
one should integrate the spectrum over the transverse 
mass and the rapidity. Hence one get for the total number 
of particle species i as 



th 



N^ = Nf" + N, 



(6) 



Starting point of particles emitted in the heavy- ion col- 
lisions is the Cooper-Frye formula[19] for the particle in- 
variant spectrum. 
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Dobler et a/. [16] have assumed the boost invariance in 
the longitudinal direction [20] and the ellipsoidal geome- 
try with azimuthal symmetry. Thus the expansion veloc- 
ity at (t, z) is vl = z/t and the corresponding rapidity 
is ry = tanh^ z/t whose value at the surface is rjmax- 
The transverse rapidity at the radial distance r from the 
center is p{r) = Po{r/rmax)°', where a = 1 is chosen 
as it gives better fit to the particle spectra. The trans- 
verse radius may depend on the longitudinal position. 



Using the above equation, one can fit the ratios among 
hadrons to get the temperature Teh, the chemical poten- 
tials pb and p^. The parameters, rimax, Po, Ro and the 
overall factor Ve/f do net affect the ratios as expected 
and their values at the chemical freeze-out cannot be de- 
termined. However, po and r]max play crucial roles in 
determining the slope of the transverse mass spectra and 
width of the rapidity distribution of the total charged 
hadrons, correspondingly. It should be noted that the 
ranges of the rapidity integration in Eq. (5), which may 
differ from particle species to another, should match the 
experimental cut-off. Due to the lack of detailed informa- 
tion we have integrated from -3 to +3 for all the particles 
and this may be one source of errors in this calculation. 
Another source of error is the lower limit of integration 
over the transverse momentum. Wc have chosen it to be 
zero but it is not. In the future studies it should match 
the experimental value. 

Depending on the data processing procedure whether 
the weak decay contribution from resonances is included 



or not, weak decay contribution should be correspond- 
ingly added or not in calculating the resonance contri- 
bution. For example, the PHENIX data for protons and 
anti-protons were already subtracted by the weak decay 
contribution from the lambda and anti-lambda particles 
whereas for pions, they arc included. Thus wc have fit- 
ted PHENIX data separately from the STAR data. In 
fitting the PHENIX ratios, we subtracted the pions from 
the decay of lamda's as subtracting out the contribution 
to protons and including it to pions is very cumbersome 
and the contribution to pions is negligible. 

The fitted parameters are listed in the TABLE 1. The 
first low is the fitted parameters from STAR data and 
the second one is for PHENIX data. The difference in 
the temperatures is less than 1 MeV while chemical po- 
tentials differ slightly. The difference may come from 
the experimental cuts for rapidity and transverse mo- 
mentum, especially the lower limit of pt, for each par- 
ticle species. The chemical freeze-out temperature of 
173Mev and chemical potentials are in agreement with 
the values from other statistical models(160-170 MeV). 
The columns which is empty in this table are irrelevant 
parameters as mentioned earlier. In the STAR analysis 
we have omitted S/tt ratio since it cannot be fitted well 
in this calculation. The value for x^/^ in TABLE 1 is 
obtained without S/tt ratio while when we include it the 
value for x^/n increases to 2.85. The (j>/K~ ratio by 
PHENIX collaboration cannot fitted in this calculation 
but the value is even twice smaller compared to STAR 
measurement, x^/^ value in TABLE 1 for PHENIX data 
is obtained without this ratio. In Fig. 1 and Fig. 2, the 
results of our chemical analysis are shown and compared 
with the STAR data and PHENIX data, respectively. 
Fig. 1 shows that the fitting is very good except for the 
S/tt ratio and in Fig. 2 fitting is extremely good for all 
the ratios except for (j>/K ratio, which differs from STAR 
data. Our fitting of this ratio to STAR data is rather 
good. 

Below the chemical freeze-out, the number of thermal 
hadrons of any type i is now fixed and known numeri- 
cally as the number of inelastic collisions is not sufficient 
enough. Whereas elastic collisions are still abundant 
to make the local thermal equilibrium and thus during 
the continuing expansion of the system, the tempera- 
ture drops until the number of elastic collisions is not 
enough to maintain even the thermal equilibrium. When 
the system breaks up the hadrons have the momentum 
distribution at this point, the so called thermal freeze- 
out. Now the transverse mass spectra and the rapidity 
distribution of measured hadrons are given from Eq. (3) 
by integrating over either the rapidity or the transverse 
momentum. The resonance contributions [16] should be 
correctly added. Explicitly, the transverse momentum 
spectrum is given as 
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and the rapidity distribution is given as 
d'^Nf^ 



dy 



mTdniTdy 



)mTdmT + (res. contr. 



(8) 



It should be stressed again that Eq. (4) does not hold 
in above equations since the chemical equilibrium is al- 
ready broken. Rather, jii can be found from the following 
equation such that it gives the right number of the ther- 
mal particles, Nf^ , which is known from the chemical 
freeze-out. 



Nj^ ^ 



rriTdrriTdy 



d^N. 



th 



rriTdmTdy 



{T, iJi, ri„iax,Po, Ro) 



(9) 

Instead of using the above equation for the number to 
calculate /i^, one can determine it relative to the pion 
chemical potential from the fixed ratios relative to pions, 
R^^ = Nl'^/N*J\ Then ii, can be obtained from 



Mi = Mtt + 7" In [Ri-, 



J J raTdmTdy(d^N^ /mTdmTdy) 



J J 'mTdraTdy{d'^N'^ /mTdmTdy) 

(10) 
where the ' denotes that exp (fXi/T) is absent in this equa- 
tion compared to Eq. (3). Now Vejf ■ e^"'^ acts as an 
overall constant for all the particle spectra. 

The resulting thermal freeze-out parameters are listed 
in the TABLE 1 together with the chemical freeze-out 
values. Thus obtained transverse momentum spectra of 
the charged pions, kaons and protons are shown together 
with the experimental data by the STAR collaborations 
in Fig. 3. The quality of fitting is rather good show- 
ing that the blast-wave model works very well for the 
transverse mass spectra in heavy-ion collisions. 

The parameter rjmax which is the rapidity of the fire- 
ball at the longitudinal surface affects very little on the 
transverse momentum spectra. However, it plays crucial 
role in determining the width of rapidity distribution and 
thus we have fitted rjmax to the pseudo-rapidity distribu- 
tion of the charged hadrons measured by the PHOBOS 
collaboration by fixing all other parameters. The fitted 
value is 5.0 as listed in the TABLE 1. In Fig. 4, we 
have drawn the rapidity distribution of charged hadrons 
together with the PHOBOS data. The shape especially 
the width of the rapidity distribution fit well to the data 
but the magnitude in the midrapidity region overpredicts 
the data. This discrepancy may come from the low pT 
cut-off in experiments while we have integrated over pT 
from zero. The wiggle near y ~ region is due to the 
errors in calculating the resonance contribution and will 
vanish if we make y-bin in our calculation smaller. It 
should be pointed that the statistical model cannot fit 
the large width of the in distribution and the present 
model is very successful in fitting the various aspects of 
the hadron production except for the correlations. As 
mentioned earlier Ve/f ■ e^^'"^ is an overall constant now. 
For Veff we get 15597, and ^^ = 115 MeV. 



TABLE I: Fitted values for each parameter: The firs and 
second lows summarize the chemical freeze-out values of the 
data by STAR and PHENIX collaborations. Empty columns 
are irrelevant parameters. 





T 
(MeV) 


fJ-B 

(MeV) 


fJ.s 

(MeV) 


Is 


f]max 


Po 


x'/n 


STAR 


173.4 


18.5 


7.9 


0.986 






1.4 


PHENIX 


173.9 


26.4 


6.0 


1.01 






0.12 


thermal f.o. 


121.1 






0.986 


5.0 


1.03 


5.3 



III. 



SUMMARY AND CONCLUSION 



in Au+Au collisions at 200 A-GeV. This is a consistent 
way of describing all the features of hadron production in 
heavy-ion collisions except for correlations. The result of 
fitting is very good except for S/tt ratio from STAR col- 
laboration and (t>/pi ratio from PHENIX collaboration. 
The parameters thus obtained are in agreement with the 
published values obtained separately either by the chemi- 
cal analysis or the thermal analysis. It should be stressed 
that in order to analyze the particle production detailed 
information on the experimental cuts for pj- and y for 
each particles is needed. Further studies with LHC data 
is underway. 



Assuming that the chemical freeze-out occurs earlier at 
higher temperature and thermal freeze-out occurs later 
at lower temperature, a blast-wave model with resonance 
contribution have been successfully applied to fit the ra- 
tios among hadrons, the transverse momentum spectra of 
pions, kaons and protons by STAR and PHENIX collab- 
orations and the rapidity distribution of all the charged 
hadrons by PHOBOS collaborations measured at RHIC 
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FIG. 4: Comparison of the rapidity distribution of charged 
hadrons by PHOBOS collaboration[25] and calculation. Only 
etamax and overall constant are fitted and all other parame- 
ters are same as chemical analysis. 



